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Abstract. We describe experiments on monodisperse spherical particles in an 
annular cell geometry, vibrated from below and sheared from above. This system shows 
a freezing/melting transition such that under sufficient vibration a crystallized state 
is observed, which can be melted by sufficient shear. We characterize the hysteretic 
transition between these two states, and observe features reminiscent of both a jamming 
transition and critical phenomena. 
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1. Introduction 

Granular materials exhibit phases analogous to conventional solids, liquids, and gases, 
in spite of being athermal and dissipative [1 . Due to the dissipation, energy must be 
supplied in order to sustain a dynamical state. Shearing and vibration are two common 
means to inject energy into granular systems. Shearing a granular material can compact 
and crystallize it 0, but also melt it [3]; tapping will compact it in thin vibrated 
layers there can be coexistence of crystallized and disordered states [Sj; and highly 
vibrated granular systems become gas-like. From a large phase space of variables we 
vary only two, the shear rate and vibration amplitude, and study the interaction of the 
two energy injection mechanisms. 

Without vibration, sheared granular materials undergo a phase transition from 
solid-like to fluid-like behavior: the particles must become unjammed (which typically 
involves dilation) before they can move. We seek to understand what effects vibrations 
have on such transitions, and on the characteristics of the states on either side of the 
transition. This is particularly interesting given that granular systems are athermal, 
and one might naively expect that vibrations would play a temperature-like role. 

We perform experiments in a classic geometry, annular shear flow lEllZlElini, with 
monodisperse particles, shown schematically in Figure H Shear and vibration provide 
competing effects, with the system evolving to a crystallized state when the kinetic 
energy provided by the vibration is greater than that provided by the shear. The 
transition is hysteretic, and fluctuations in the packing fraction and the breadth of the 
force distribution both become large as the crystallized state is approached, in similarity 
to phase transitions in other systems. 

The physical parameters that characterize the system include the amplitude A and 
frequency / of vibration, the height H and mean radius r of the annular container, 
the diameter d and the density p of the particles, the rotation rate Q of the upper 
shearing surface, and the mean pressure P on the layer (here characterized at the base 
of the layer). From these physical parameters, it is possible to define other dimensioned 




Figure 1. Schematic cross-section of experiment (not to scale). 
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r vibrational acceleration / gravitational acceleration 

9 

H — cell height / particle size 

(Ji 

I — particle scale velocity / acoustic velocity 

c 

J apparatus scale velocity / vibration velocity 

Au 

K vibration velocity / acoustic velocity 

c 
j[ 

L — vibration length scale / particle length scale 
a 

M — shear time scale / vibration time scale 

UJ 

P 

N applied pressure / hydrostatic pressure 

pgd 

r 

R — cell radius / particle diameter 
a 

Qr 

f2 shear velocity / particle velocity 



gd 



Table 1. Dimcnsionless ratios, with uj = 27r/, 7 = ftr/h, and c = \/P/p 

parameters, such as the shear rate, 7 = Qr/h, as well as a number of dimemsionless 
parameters which we list in Tabled 

r, / and J are three key parameters from the list in Tabled Of the ten listed, there 
are only seven independent parameters: for example, F = K'^N/ L and the four velocity 
ratios {H, I, J, K) only represent three parameters. In the experiments described here, 
we have fixed /, P, and and therefore only explored a small region of the available 
phase space. 

2. Experiment 

The experimental apparatus consists of an annular region containing nearly 
monodisperse polypropylene spheres of diameter d = 2.29 to 2.39 mm and density 
p = 0.90 g/cm^, as shown in Figure ^ with the pressure P and volume V (height h) 
set from below by a spring within an electromagnetic shaker. The particles are sheared 
from above and vibrated from below, while the sidewalls are stationary. A more detailed 
description of the apparatus is given in To characterize the states, we obtain high- 
speed video images of particles at the outer Plexiglas wall, laser position measurements 
of the bottom plate (cell volume), and force time series from a capacitive sensor flush 
with the bottom plate. For the experiments described in this paper, we fix the frequency 
of vibration (/ = 60 Hz) and number of particles (A^ ~ 71200), and vary the amplitude 
of vibration A and shear rate Q. We vary the nondimensionalized peak acceleration 
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Figure 2 

and n = 0.078 



Sample images, viewed from outer wall, (a) Crystallized state: F — 2.0 
http : / /nile . physics .ncsu.edu/pub /cryst . mpg (b) Disordered 



Movie at 



state: F = 2.0 and Q, ~ 0.47. Linear (L) and hexagonal (H) clusters marked by black 
boxes. Movie at http: / /nile. physics. ncsu.edu/pub/disord. mpg" (c) Phase diagram for 
crystallized and disordered states as a function of Cl and F. Dashed line is J = 1. 
Adapted from 



r = A{27if y/g from to 6, and nondimensional shear rate fl = Qr/ y/gd from 0.058 to 
9.3. 

3. Description of States 

In the regime < f2 < 10 and < F < 6 we observe two distinct granular states 
of matter: crystallized and disordered. Sample images and movies of these two states 
are shown in Figure |2t'b, as viewed from the outer wall. For < 1 and F < 4, we 
observe that the phase boundary between the two states roughly corresponds to a curve 
where the characteristic velocities of the two motions are equal. This corresponds to 
the dimensionless number J of Table 
Qr 

'Va « 

and the curve J = 1 is shown by the dashed line in Figure |2t. Below, we characterize 
these two states, with further details to be found in [3]. 

Crystallized State: In the solid-like state (see Figure 12^-), the balls crystallize into 
a hexagonally close-packed configuration {i.e. a 3D crystalline structure) here visible 
only at the outer wall although the order persists across the layer. The contact between 
the upper layer of the granular material and the shearing wheel is intermittent, with 
stick-slip motion of the top ~ 2 layers in the manner of JU]- The distribution of forces 
measured at the bottom of the layer is bimodal (see Figure E)); this indicates that the 
material is responding as a solid body moving up and down with the sinusoidal vibrations 
of the bottom plate. 

Disordered State: In the disordered state, some order remains in the form of 
hexagonally-packed clusters and linear chains of particles at the outer wall, as marked 
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Figure 3. Force probability distribution functions for three values of f2 at F = 2.0. (a) 

on linear scales, dimensioned and (b) on log-linear scales, normalized by mean force. 



in Figure |2b- For states with Q well above the transition, linear chains dominate over 
hexagonal clusters, with both existing intermittently throughout the disordered regime. 
These chains may correspond to the planar ordering reported recently by Tsai et al. 
121 E]- The velocity profile extends deeper into the layer (in the vertical direction) 
than in the crystallized state. Force distributions measured at the bottom plate show 
the exponential-like tails characteristic of many granular experiments in disordered, 
unvibrated granular materials (see Figure Ej). They also fall to zero at low force, as seen 
in earlier experiments by Miller et al. [7]. 

For a geometrically similar system, but unvibrated and exposed to a compressional 
force, shear ordered the system into horizontal planes of hexagonal packing, each slipping 
past the others |21 HI] • Such a state is different from the 3D crystallized state observed 
here, in which the layers in the bulk are stationary with respect to each other. An 
interesting question is how shearing creates order or disorder depending on the presence 
or absence of vibration. A useful way to distinguish the ordered and disordered regimes 
is via the ratio 



(2) 



involving the shear time scale 'j = ilr/h to the acoustic time scale, d^pjP ^ calculated 
from the pressure P and density p jT2]. The experiments described in this paper have 
pressures of around P = 20 Pa and shear rates of 7 = 0.3 to 40 Hz, leading to / = 5 x 10^^ 
to 0.75. In j2], the shear rates are slower, 7 = 0.05 to 0.5 Hz, and pressures are higher, 
P = 2000 Pa, so that the system is clearly in the quasistatic regime with / = 3 x 10~^ 
to 3 X 10"^ 

While it is perhaps surprising that we find as simple a result as a phase transition 
at J ~ 1, the presence of these other important control parameters give hints into the 
breakdown of crystallization for large P. Figure|2t shows that for Cl ^ 0.7, crystallization 
was not observed above P = 4, possibly indicating the re-emergence of disorder due to 
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Figure 4. Characterization of velocity profiles at F = 2.0 and various fl. (a) 
Azimuthal velocity measured at outer wall as a function of depth measured. Lines 
are fits to Equation Q. (b) Solid body rotation function of f2, (c) decay 

length j/o as a function of fl and (d) slip coefficient a as a function of CI. 

granular-gas like behavior. Further experiments varying (L, M, N) will be necessary to 
discover which of these determine the high-F boundary of the crystallized phase. In 
addition, the parameter H controls finite size effects. 

4. Shear Localization 

Granular materials commonly exhibit shear banding, with an exponentially decaying 
velocity profile away from the shearing surface. As shown in Figure |3^, this shear band 
behavior is seen in both the crystallized {fl = 0.087) and disordered states (f2 = 0.87) 
for F = 2.0. To obtain the velocity profiles, we tracked individual particles visible at the 
outer wall using a high-speed video camera and determimed trajectories for each. We 
then binned the resulting velocity components by depth to construct velocity profiles. 
We characterize the azimuthal velocity, v{y), by fitting the profile to the form 

v{y) = v^ + aviie-y/y" (3) 

where Voo is the solid-body motion at the bottom plate, fy is the known azimuthal 
velocity of the shear wheel at the outer wall, a is the efficiency with which that velocity 
is transmitted to the top layer of granular material, and yo is the decay length of the 
velocity with depth. 
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Figure 5. Volume V of cell, scaled by minimum observed volume Knm as a function 
of (a) ri (at r = 2.0) and (b) F (at 17 = 0.27). Triangles point in direction of steps. 



In the crystallized state, the shear is localized almost entirely to the first layer of 
particles (small i/o) , while in the disordered state the shear band extends several particles 
into the layer. The slip at the upper plate is lowest in the disordered state, where the 
uppermost particles are in constant contact with the shearing wheel. Note that the 
system is more dilated in the disordered state, and has a larger pressure j^. While the 
disordered state has greater slip (foo) at the bottom plate, the scaled slip values {voo/v\\) 
are in fact lower than in the crystallized state, visible in Figure|3^. For disordered states 
with Q > 0.4, the shear bands appear to have reached a steady state since they are all 
parameterized by the same values. 

5. Transition 

We examined the transition from the disordered to the crystallized state by first 
preparing a disordered state at high Q, then adjusting Q to the value of interest. We 
then performed two runs, one at constant Q = 0.27 (starting from F = 0) and the other 
at constant F = 2 (starting from Q = 8.4). The mean volume measured for each step of 
these two runs is shown in Figure 

For steps of decreasing Q (Figure Et^) the system compacts logarithmically until 
reaching Qc, after which the system undergoes a first order phase transition to the 
crystallized state. Below only a small amount of additional compaction occurs, 
to a state with a volume Vmm, for which the packing fraction is = 0.69. When 
Q is increased, the transition back to the disordered state is hysteretic, occurring for 
^Ih ^ 2fic- 

For steps of increasing F (Figure |S|d) the system also compacts. However, runs 
approaching the transition are difficult to repeat quantitatively, since there is a great 
deal of intermittency in the cell volume (see Figures Eb, and jH] for details). For 
F > Fc the system is in the crystallized state. The transition also appears to be first 
order, but in this case the hysteresis is so extreme that the material was not observed 
to re-expand when we decreased F. 

We wish to understand why a crystallized state can disorder by increasing fl, but 
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Figure 6. Characteristics of force probability distributions as a function of fl at 
r = 2.0. Triangles point in direction of steps in fl; solid points are crystallized phase; 
dashed line is guide to the eye. 



not by decreasing T. In the case of increasing the stick-shp behavior in the top layers 
of the crystalhzed state is affected by the speed of the upper shearing wheel. As Q 
increases, more horizontal momentum is transfered to the upper layer of balls, which 
results in longer regions of flowing particles. Eventually, the whole layer can be seen 
to mobilize and the disordering begins to take place throughout the cell. In contrast, 
for increases in A (and hence F) no such increased momentum transfer takes place, 
and the results are similar to the irreversibility observed for compaction by tapping |3]. 
This transition shows some similarity similarity to the "freezing-by-heating" transition 
seen in in which individual particles with tunable noise are seen to crystallize as 
their noise level is increased. Such a system also shows hysteresis in returning to the 
disordered, mobihzed state. 

For the run at F = 2.0, seen in Figure EK with steps downward in Q, we observe 
signatures of the phase transition from disorder to crystallization via both the volume 
fluctuations and the force distribution, as shown in Figures IHl and [3 As Q ^ Qc 
from above, both the volume fluctuations (measured from the variance of V{t)) and the 
breadth of the force distribution (measured by the kurtosis, or fourth scaled moment, 
of F{t) on the force sensor) become large. In addition, the first order nature of the 
transition is visible in other characteristics of the force distribution, such as the mean, 
standard deviation, and skewness (see also Figure Ej). 

In granular systems, Edwards and coworkers J3] have introduced a temperature-like 
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Figure 7. Volume fluctuations as a function of ft, scaled by average volume Vq. Inset: 
Volume fluctuations as a function of Vo/Knm at each 17. Triangles point in direction 
of steps in 17; solid points are crystallized phase; dashed line is guide to the eye. 



measure, the compactivity, defined as X = {dV/dS)N by analogy with thermodynamics. 
The central idea is that lower packing fractions correspond to a greater freedom for 
particle rearrangement, and hence a higher compactivity. In the statistical mechanics 
of ordinary second order phase transitions, susceptibilities such as {d'^A/dT^)v (for free 
energy A) are singular at the critical point. For example, the specific heat at constant 
volume is Cy = {dE/dT)v = —T{d'^A/dT'^)v- When described in terms of fiuctuation- 
dissipation relations, ksT'^Cv = {{E — EoY), where E is the energy of the system and 
£^0 its mean value. One expects energy fluctuations, and hence Cy, to be singular at the 
critical temperature Tc. By contrast, at a first order transition, discontinuities occur 
in densities, but one does not expect divergent fiuctuations. Since V has taken the 
place of E in the Edwards formalism, the hallmark of a critical transition is increased 
fiuctuations in the volume of the system as we approach Xc, the critical compactivity. 
In our experiments, volume (and hence X) is set by Q, and the inset to Figure [7| shows 
apparently singular behavior for the volume fiuctuations as a function of the volume. 
It is interesting that in these experiments, we see a discontinuity in the density, but 
also an indication of a singularity in the volume fiuctuations. The magnitude of the 
fiuctuations observed in the disordered state is similar to those in observed in jl5j . 
where the standard deviation of the packing fraction is approximately 10~^. 



6. Intermittency 

The apparently singular volume fiuctuations near J = 1 come from the fact that the 
system exhibits intermittency in its state. The system is in fact spatially inhomogeneous, 
with instances of small V being crystallized in the majority of the cell and instances 
of large V being majority disordered. By examining the properties of the system in 
this intermittent regime, we are able to compare a broad range of states for nearly the 
same parameter values. The only varying parameters are the volume and pressure of 
the system, which are related to each other by a proportionality constant due to the 
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Figure 8. Characterization of velocity profiles in intermittent regime at f2 = 0.27 
and r — l.O. (a) Azimuthal velocity measured at outer wall as a function of depth 
measured for compact (Q, t = 3.8 hrs) and dilated (□, t — 4.7 hrs) states. Lines are 
fits to Equation (PJ. (b) SoHd body rotation Vqq BiS 8j function of cell height, (c) decay 
length yo as a function of cell height and (d) slip coefficient a as a function of cell 
height. Solid symbols correspond to data from (a). 



spring constant of the shaker T\. 

The inset to Figure shows a time series of the cell height for such a run, with 
r = 1.0 and fl = 0.27. The system started from a dilated state and progressed to 
a majority crystallized state before re-dilating and re-compacting to an even more 
crystallized state over the course of approximately 10 hours. We again obtain velocity 
profiles at the outer wall, in this case while simultaneously monitoring the position of 
the bottom plate. 

Figure IHK shows two velocity profiles from compact and dilated states. Because the 
system is spatially inhomogeneous, the particles in view of the camera may in fact be 
either disordered or crystallized at any given time during these measurements, regardless 
of the height of the cell. Importantly, the fit parameters in Figure IHb-d show the same 
trends as those in Figure 0] when rotation rate fl is taken as a proxy for cell height. In 
both cases, we observe a continuum of states as the system moves between crystallization 
(compaction) and disorder (dilation). Again, it is interesting to note that the volume and 
pressure fluctuations are associated with the formation and melting of ordered clusters 
over time. Such behavior is characteristic of near-critical behavior. By contrast, at a 
thermodynamic first order transition, we would not expect to see persistent dominant 
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fluctuations. 
7. Discussion 

The two characterizations of a transition in the system we discuss above provide 
contrasting, but complimentary, information about the nature of the crystallizing phase 
transition in sheared and vibrated granular materials. The canonical hallmark of a 
transition to a jammed/glassy state is the continuous growth of the viscosity. Glass 
transitions do not in general contain first-order-like signatures, such as discontinuities in 
the volume or specific heat jTHI. For sheared colloids, there are large stress fluctuations 
near a jamming transition jT2j, and in simulations of Lennard- Jones particles, force 
PDFs are observed to broaden ^SJ ^j. Similar behavior is observed in this system 
as well, but with a density discontinuity. By contrast, jammed/glassy states are all 
disordered, while the granular system described in this paper makes a transition to a 
crystallized state. In both the glassy and crystallized cases, however, the final states are 
unable to rearrange. 

We observe similarities to critical phenomena in the increased volume fluctuations 
near the transition, a hallmark at odds with a glass transition. These fluctuations 
are similar to the density fluctuations observed at the liquid-gas critical point, which 
occur at diverging length scales. Therefore, further investigations into the nature of 
this transition should examine what length scales and order parameters are present, 
including a determination of the sizes of clusters and the spatial correlations between 
forces. Finally, we have introduced a number of dimensionless control parameters whose 
effects remain to be investigated in future studies. 
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